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Abstract 

We describe the simulation of tau decays in the Herwig++ event generator, which 
includes sophisticated modelling of the hadronic currents and full treatment of spin 
correlation effects. The structure of the simulation makes it easy to add new models 
of tau decay, change the parameters of the existing models, and use the models from 
tau decay for the decay of other particles. The results are compared in detail with an 
existing simulation, and the benefits of the new structure are illustrated by consid- 
ering the decay Xi ^ Xi + hadrons in Anomaly-Mediated SUSY-Breaking (AMSB) 
models. 
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1 Introduction 

The use of Monte Carlo event generators is an essential part of all experi- 
mental analyses, both in interpreting data from existing experiments and in 
the design and planning of future experiments. The crucial role that Monte 
Carlo simulations play in experimental studies mean it is imperative these 
simulations are as accurate as possible. 
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While the existing Monte Carlo event generators have been highly successful 
over the last twenty years, a new generation of programs is necessary for the 
LHC The reasons for this are twofold: a number of new ideas to improve 
the accuracy of the simulations have been suggested, e.g. P^PfSPfS] . and the 
existing code structures required major redesign to allow new theoretical de- 
velopments to be incorporated and to make long-term maintenance easier. 
Given the changing nature of computing in high energy physics, the natural 
choice is to write these new programs in C++. In preparation for the LHC 
a major effort is therefore underway to produce new versions of established 
simulations [HIT], as well as completely new event generators [8j in C++. 

As part of the process of writing the new Herwig++ event generator |7| we 
wish to improve many aspects of the simulation process. One area where major 
improvements are needed is in the simulation of tau lepton decays. 

In FORTRAN HERWIG the r decays were treated in the same way as the decays 
of the hadrons. The weak V — A matrix element was used for the leptonic 
decays, and the other decays were generated using a phase-space distribution 
for the decay products. In addition, there were two interfaces to the specialised 
TAUOLA [9][T0] package for r decays. The first interface [1], which was part of 
HERWIG, was capable of generating longitudinal correlations in r decays for all 
taus produced in the perturbative part of the event, while the second external 
interface [LQl was capable of generating the longitudinal spin correlations of 
taus produced in W, Z/j* and if^ decays and the full correlation effects in 
neutral Higgs decays. 

In Herwig++ we wanted to: 

(1) include the matrix elements for r decay as an integral part of the sim- 
ulation, in order to both have a unified treatment of all decays and to 
remove the inherent problems in interfacing to external packagesPI 

(2) include the absent transverse spin correlations for all r decays; 

(3) make it easy to change parameters in the hadronic currents used in the 
tau decays and the modelling of individual decay modes. 

In this paper we will describe the simulation of r decays in Herwig++. The aim 
of the simulation is to give a good description of tau decays with all the ex- 
perimentally observed decay modes with branching ratio above the 5 per mille 
level included together with a reasonable model of the kinematics of the de- 
cay. Where possible we have used models which have been compared with, or 
tuned to, experimental data. In all cases the parameters of the models can 
be easily adjusted so that they can be tuned to new experimental data as it 
becomes available. 

^ In the FORTRAN simulation there were problems with the HERWIG interface due 
to changes in the TAU OLA program. 
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First we describe the factorization of the matrix element for r decays and 
the structure of the code. This is followed by a description of the different 
models of the hadronic current which are used for the decay modes, together 
with comparisons with previous results. We then discuss our choice of models 
for the various r decay modes. One major advantage of the new structure is 
that the hadronic currents can be used for applications other than r decays 
which we illustrate by considering the weak decay of charginos to neutralinos 
in Anomaly-Mediated-SUSY Breaking (AMSB) models, where there is a small 
mass difference between the neutralino and chargino. Finally we present our 
conclusions. 



2 Formalism 

The matrix element for the decay of the r lepton can be written as 
Gf 

M = -^L^J'', L^ = m(p^J7^(1-75)m(p,), (1) 

where Pr is the momentum of the r and p^^ is the momentum of the neutrino 
produced in the decay. The information on the decay products of the virtual 
W boson is contained in the hadronic current, J^. This factorization allows us 
to implement the leptonic current for the decaying tau and the hadronic 
current separately and then combine them to calculate the r decay matrix 
element. 

In Herwig++ the most important part of the simulation of hadronic decays is 
handled by the Decayer class. This class is responsible for: 

(1) generating the kinematics of the decay products; 

(2) inserting the spin-unaveraged matrix element used to perform the de- 
cay and the wavefunctioniiJ for the decay products into the Herwig++ 
structure. 

All of the r decays are handled by the Tau Decayer class. This class inherits from 
the Decay Integrator class. The Decay Integrator includes a sophisticated multi- 
channel integrator to perform the phase-space integration of the decay modes, 
leaving the implementation of the matrix element calculation and definition 
of the resonance structure of a multi-body decay mode the only tasks which 
need to be performed when implementing a new decay model. 

The hadronic currents for a large range of modes are implemented as described 

^ In Herwig++ these are the spinors and polarization vectors of the fermions and 
vector mesons. 



3 




b) 6' 





Hw++ H 

Hw++ A ^r'' 

/ 








— 





















3,000 3,025 3,050 3,075 3,100 3,125 



Fig. 1. Correlations between the decay products for {H,A} t^t^ followed by 
r TTVr for a Higgs boson mass of 120 GeV. This plot shows the variables considered 
in [TIJ. Figure (a) shows the angle between the planes of the two tau decays and 
(b) shows the angle between the pions in the rest frame of the Higgs boson. 

in Section [3l The classes implementing these currents all inherit from WeakDe- 
cayCurrent. The TauDecayer combines the hadronic current with its internal 
calculation of the leptonic current to compute the matrix element, and gen- 
erates the momenta of the decay products using the features of its parent 
Decay Integrator class for the generation of the phase space. 

Given the spin-unaveraged matrix elements, the algorithm described in pTfT^fT 
is used in Herwig++ to include spin correlation effects in all stages of the event 
generation process. The Herwig++ structure, in particular the DecayHandler 
class, is responsible for passing the correct information between the Decayer 
objects in order to generate the spin correlations. 

Examples of these correlations for the decay of scalar and pseudoscalar Higgs 
bosons to r+r~ followed by the decays of r — > ttv^- and r — > pVj. are shown in 
Figures [1] and [2] respectively!!] In the absence of the correct correlations the 
0* distributions in both cases would be flat. The approach of p!T|12lll3|[T] is in 
good agreement with the results of TAUOLA where the two methods can be 
comparec0 but can generate the full correlations regardless of how the tau is 
produced. 



^ In all the plots version 2.1 of Herwig++ is compared with the TAUOLA-PHOTOS 
package of October 2005. 

^ The full effects can only be compared in neutral Higgs decay as this is the only 
case in which TAUOLA implements the full correlations. However, there are many 
cases where the transverse correlations are not important and for which the two 
approaches are in good agreement, see for example Section 4.4 of [15] . 
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Fig. 2. Correlations between the decay products for {H, A} t^t^ followed by 
r pvr for a Higgs boson mass of 120 GeV. This plot shows the variables considered 
in [16]. The angle between the decay planes of the two p mesons in the p'^ p~ rest 
frame is shown for (a) yiy2 > and (b) yiy2 < 0. The variables yi^2 are defined in 
terms of the energies of the particles in the respective r rest frames, yi = 

for the decay and ?/2 = a-e" P~ decay. 



In Figures [T] and [2], and all subsequent, figures the TAUOLA result is not shown 
as it is virtually indistinguishable from the Herwig++ result due to the high 
statistics used. Similarly where present the lower panel shows the difference 
between the TAUOLA [9] and Herwig++ results in terms of the statistical error 
for mass distributions, and the fractional difference for running widths. 



3 Hadronic Currents 



We have implemented a wide range of models for the weak hadronic currents. 
The various models are described in this section together with comparisons 
with the results of TAUOLA, where appropriate. 



3. 1 Pseudoscalar Meson 



The simplest hadronic current is that for the production of a pseudoscalar 
meson, e.g. the current for the production of vr^ in the decay of the tau. The 
hadronic current can be written as 



J' = fpp'p, 
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Mode 


Herwig++ 


TAUOLA 


Difference 




rpartial/10-13 GeV 


rpartiai/10-13 GeV 


/lO-i^GeV 


PseudoScalar Meson 


vr 


2.4334 


2.4334 





K- 


0.16611 


0.16611 





Two PseudoScalar Mesons 


via Intermediate Vector Mesons 




5.3998 ± 0.0002 


5.3998 ± 0.0002 


0±3 




0.081295 ± 0.000002 


0.081293 ± 0.000001 


0.02 ± 0.03 




0.156196 ±0.000004 


0.156195 ±0.000003 


0.01 ±0.04 




0.0024248 ± 0.000001 


0.0024251 ± 0.000001 


-0.03 ± 0.02 


Ktt via Intermediate Scalar and Vector Mesons 




0.081292 ± 0.000003 


0.081293 ± 0.000001 


-0.01 ± 0.04 




0.156200 ±0.000005 


0.156195 ±0.000003 


0.05 ±0.06 


Charged Lepton and Neutrino 


e I'e 


4.0491 ± 0.0002 


4.0492 ± 0.0002 


-1±3 




3.9380 ± 0.0002 


3.9380 ± 0.0002 


0±3 



Table 1 

Partial widths for two- and three-body decay modes of the r calculated with Her- 
wig±± and TAUOLA. In order to compare the results of Herwig±+ with TAUOLA, 
the masses of the decay products in TAUOLA were adjusted to the Herwig±± val- 
ues, as were the pion and kaon decay constants and the Cabibbo angle. For the two 
meson, via intermediate vector mesons, decays the parameters of the resonances in 
the form-factors in Herwig±± were set to those used in TAUOLA, the CPC version 
was used for the tt'tt^ and K~ modes and the CLEO version for the K-k modes. 
For the Ktt modes including scalar resonances a transverse form of the projection 
operator was used together with the resonance parameters from the CLEO version of 
TAUOLA in Herwig±±. In the leptonic decays the radiative corrections in TAUOLA 
were switched off. For the Kt: and K~ modes our choice of the normalisation 
was used. 



where pp is the momentum of the pseudoscalar meson and fp is the pseu- 
doscalar meson decay constant. 

The partial widths for the decays r~ vr"//,- and K~Vt- are given 

in Tabled] and agree with those from TAUOLA. Another important check is 
that the spin correlation effects are correctly implemented. The correlations 
in Higgs decays followed by the decay of the tau to a single pion are shown in 
Fig. [Hand are in excellent agreement with the results of TAUOLA. 
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Fig. 3. The differential decay rate with respect to the mass of the lepton- neutrino 
pair for t~ l~i>iVr. The TAUOLA result was generated with the radiative correc- 
tions switched off. 

3.2 Vector Meson 



The current for the production of a vector meson is given by 

= V2gye*yr (3) 

where ey is the polarization vector for the outgoing meson and gy is the decay 
constant of the vector meson. 

This current was included to test some aspects of the spin correlations, test the 
treatment of off-shell effects and for future extensions, but is not used in tau 
decays because all such decays are better modeled by a complete description 
of the the production and decay of the vector meson as implemented in the 
various models described below. 



3.3 Charged Lepton and Neutrino 



The current for weak decay to a lepton and the associated anti-neutrino is 
given by 

J^ = n(p,)7ni-75)^(pp), (4) 
where is the momentum of the anti-neutrino and pi is the momentum of 
the charged lepton. 

The mass distribution of the C.~v^ produced in the decay — s> I^viVt is shown 
in Fig.O As can be seen there is good agreement between Herwig++ and 
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Fig. 4. The mass spectrum for tt^tt^ in the decay r -k^-k^Ut- a) shows m^_^o 
compared with CLEO data taken from [19] and b) shows m^-^o compared with 
Belle data from [20j. The results of Herwig++ using the form-factor parameters 
from the fits of CLEO LlSj, Belle [20j and ALEPH [21j are shown together with the 
result from Herwig++ using the TAUOLA parameters. 

TAUOLA in terms of both the differential distributions and the partial widths 
given in Tabled! For these comparisons the masses of the particles in TAUOLA 
were set to the Herwig++ values and the radiative corrections switched off. 

In general we have not yet considered radiative corrections to tau decays in 
Herwig++. However, we do have a treatment of electromagnetic radiation in 
1^2 decays [17j based on the YFS formalism [18] which resums the dominant 
soft radiation to all orders and treats the universal large coUinear logarithms. 
This can be applied to r decays in the approximation that the decay is treated 
as a series of 1 ^ 2 processes, which is reasonable for many of the decay 
modes. The YFS formalism can be systematically improved to incorporate 
exact, process specific, higher order corrections so in the future we can extend 
this simulation to include the full radiative corrections in tau decay. 



3.4 Two Pseudoscalar Mesons via Intermediate Vector Mesons 



The weak current for the production of two mesons via the p or K* resonances 
has the form 

J' = (pi - P2)u (g"" - ^] ^ E c^kB.iq'), (5) 
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Fig. 5. The mass spectrum of a) K'tt^ and b) -k^K^. In a) the solid line is the 
result of the model described in Section [3.41 using the TAUOLA parameters for the 
form- factor and in b) the solid line is the result of the model described in Section [331 
with no scalar component, using with a transverse form of the projection operator 
and the parameters of the vector form factor set to the TAUOLA values. In both 
plots the dashed line was generated using the model from Section 13.51 using the fit 
of 



where pi,2 are the momenta of the outgoing mesons, q = Pi +P2, Bk{q^) is the 
Breit-Wigner distribution for the intermediate vector meson k and ak is the 
weight for the resonance. The Breit-Wigner terms are summed over the p or 
K* resonances that can contribute to a given decay mode. 

The models of either Kiihn and Santamaria [22], which uses a Breit-Wigner 
distribution with a p-wave running width Eqn. lA.2] or Gounaris and Saku- 
rai ^23j, which uses the form given in Eqn. lA.4l are supported for the shape of 
the Breit-Wigner distribution. 

The following decay modes: —>■ tt^tt^; K*^ —>■ K~7C^; K*~ K^n^; 
p~ — >• K~K^; K*~ K~r] are supported with weights, W\ \/2; 1; 1; 

^. respectively!!] 

There has been a number of recent studies of the tt^tt" mode p^21|20j . The fits 
from these studies are compared with the data from CLEO [12] and Belle [20] 
in Fig. m together with a comparison of the results of Herwig++ and TAUOLA. 
The partial width using the TAUOLA parameters is compared with that from 
TAUOLA in Tabled The CLEO, ALEPH and Belle fits differ in the amount 



^ It should be noted that this normalisation for the Kt: and KK modes is different 
from that in TAUOLA |9|lUj . However, it agrees with that of [25J for the Kt: modes 
and [26] for the Krj mode. 
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Fig. 6. The mass spectrum for a) K" and b) K^rj. For the K~ spectrum the 
sohd hues were generated using the parameters from the CPC version of TAUOLA 
for the form factor while the dashed hne uses the default Herwig++ parameters for 
the rho form factor. 

of the p(1700) resonance which is present, giving destructive interference in 
the high mass region. We have therefore chosen to use the parameters of the 
CLEO fit using the model of Kiihn and Santamaria [22j, which lies between 
the ALEPH and Belle fits, as our default choice. 



There have been a number of recent studies of the Kti modes P7ll28f24j which 
favour a low-mass enhancement and little contribution of the ii'*(1410). We 
have therefore chosen use a model which includes scalar resonances, as de- 
scribed in Section 13. 5[ by default for the Kti modes. For this model we use 
same parameters as the CLEO version of TAUOLA, which includes the i^'*(892) 
and i^'*(1680) resonances. The mass spectrum for the K^n^ mode shown in 
Fig. [5^ shows good agreement between Herwig++ and TAUOLA, as does the 
partial width for the Ktt modes given in Table [H 

The mass spectrum for K~K^ generated by Herwig-|-+ is compared with that 
generated by TAUOLA in Fig. [6k, again there is good agreement between Her- 
wig+-|- and TAUOLA for both the shape of the distribution and the partial 
width for the decay mode given in Table [H The spectrum for the Krj mode is 
shown in Fig.[6]D. 



3.5 Ktt via Intermediate Scalar and Vector Mesons 



Unlike the vr+vr^ decay of the tau the Kir decay mode can occur via either 
intermediate scalar or vector mesons. We therefore include a model for the 
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0.5 1.0 1.5 2.0 2.5 3.0 

Q/GeV^ 

Fig. 7. The hadronic structure function Wtot for r — > K'kVt-. The solid black line 
shows the result with the default parameters of [29] , the short dashed line shows the 
result without the scalar contribution, the dot-dashed line shows the result without 
the i^*(1410) and the dotted hne shows the result with neither the i^*(1410) or 
scalar contribution. This plot is the equivalent of Fig. la of |29j . 

current for the Ktt decay mode including the contribution of both vector and 
scalar resonances based on the model of [29j. The current is given by 

(6) 

where pi^2 are the momenta of the outgoing mesons, q = Pi + P2, BWfc(g^) is 
the Breit-Wigner distribution for the intermediate mesons and is the weight 
for the resonance. The sum over the resonances is over the vector K* states 
in the first, vector, part of the current and the excited scalar K* resonances 
in the second, scalar, part of the current. By default the vector part of of the 
current includes the i^'*(892) and i^'*(1410) states and the scalar part of the 
current includes the Kq{1A30) together with the option of including the fi;(800) 
to model any low-mass enhancement in the mass of the Kir system, although 
additional resonances can be included if necessary. 

The mass spectrum for the n'K^ mode is shown in Fig.[5b, where the pa- 
rameters have been chosen to give the same form factor as TAUOLA. There is 
good agreement with TAUOLA for both the shape of the distribution and the 
partial widths, which are given in Table [H 

The results for the hadronic structure function defined in [29] are shown in 
Fig. [7] for different values of the parameters. This shows good agreement with 
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Fig. la of [29], upon which it is based. The branching ratios for the different 

parameters values are also in good agreement with the results of |29| . 




The most recent study [21] of the Kn mass spectrum sees a low-mass enhance- 
ment which is modelled using the scalar k resonance. By default we use a set 
of parameters this model which reproduces the fit of [23] to the mass spectrum 
with the K, i^*(892) and i^'*(1410) resonances, as shown in Fig.[5]for both the 
K~TT^ and tx~K^ mass distributions. The low- mass enhancement due to the k 
can clearly be seen in these distributions. 

3.6 Three Pseudoscalar Mesons 

In order to simply the implementation of a number of standard currents for 
the production of three pseudoscalar mesons we define the current in terms of 
a number of form factors. The current is defined to be [9] 



J'=[9"'~^-^] - P,Y + F2{p, - PiY + Fz{Pi - P2Y] (7) 



where pi,2,3 are the momenta of the mesons in the order given below and 
-Fi_^5 are the form factors. We use this approach for a number of three me- 
son modes which occur in r decays: tt~tt~tt^\ tt^tt^tt"; K~tt~K^\ K^tt~K^\ 



K-TT^K^; TT^Ti^K-; K~'r-tx+- ti-R'^ti^; n-n°r]; K'^s^-R^s', KI-k-RI; R^g^-R^. 



The current is implemented in terms of these form factors in a base class so 
that any model for these currents can be implemented by inheriting from this 
class and specifying the form factors. 

We currently implement three models for these decays, the general model 
of [22]l30f9] which treats all the decay modes, the model of CLEO [31] for the 
three-pion modes and the model of [25] for the kaon modes. 

3. 6. 1 General Model 

This is the implementation of the model [22][30ll9] which uses the form of [22] 
for the Oi width. The form factors for the different modes are given in [30 .9]. 
A Breit-Wigner distribution, Eqn. lA.2] with a running width is used for the 
ai. As the dominant decay of ai is to three pions, a more complicated form 
of the running width is used. A parameterisation of the energy dependence of 
the running width is given [22] for the default parameter values for the p form 
factor. Instead of using this parameterisation we calculate an interpolation 
table for the running width at initialisation with the actual p form-factor 





+g^F4 + iF5e'""^^py^pl 
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Fig. 8. The running ai width calculated using the model of ^22j. a) shows the running 
width for the parameters used in the Section 13.6.11 and b) shows the running width 
for the parameters used in the Section 13.6.31 The result from numerically evaluating 
the matrix element in Herwig++ is shown. 

parameters. Our calculation of the running width is compared with the result 
of [22] in Fig. [HI there is good agreement between the two results. 

The results for the partial widths for all the modes are compared with the 
results from TAUOLA in Table [2] which shows good agreement between the 
two programs. The mass distributions for the hadronic system and the tt+tt^ 
subsystem, which contains the p resonance, for the r~ 7r~7r"7r+z/T- mode are 
shown in Fig.O The mass distributions of the tt^tt^tt" system and the 7r~7r° 
system are shown in Fig.[TUl for the r~ —>■ 7r°7r°7r^t^^ mode. 

The mass distributions of the K~tt~K^ system and the 'k~K^ system are 
shown in Fig.[TT]for the r~ — > K~tt" K^Ur decay mode. The mass distributions 
of the K^n'K^ hadronic system and the 7t~K^ system are shown in Fig.[T2] 
for the T" K^TT~ K^Vj. decay mode. The mass distributions of the hadronic 
system and the K~tt^ subsystem are shown in Fig.[T3]for the r~ —>■ K^tt^K^Ut 
decay mode. 

The mass distributions of the hadronic system and the 7i^K~ subsystem are 
shown in Fig.[T4]for the decay mode r~ Tc^ir^K^Ur- The mass distributions 
of the K^Ti^Ti^ system and vr+vr" subsystem are shown in Fig. [15] for the 
decay K^Ti^-K^Vr- The mass distributions of the vt^A'^tt^ system and 

the 71^ subsystem are shown in Fig. [TBI for the decay -k^ K^ti^Vt-. 

Finally the mass distributions of the ti'tt^t] system and the h'tt^ subsystem 
are shown in Fig.[T7|for the TX^TX^rjVj. decay. It should be noted that this 

mode is very sensitive to the value of the rj mass used. 

In all cases there is good agreement between Herwig++ and TAUOLA. 
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Mode 


Herwig++ 




TAUOLA 




Difference 




rpartial/10-13 GeV 




rpartiai/10-13 GeV 




/lO-i^ GeV 


General Model 


vr vr vr^ 


1.4632 ± 0.0001 




1.4633 ± 0.0001 




-1 ± 1 


n n 
vr vr^TT^ 


1.4957 ± 0.0001 




1.4958 ± 0.0001 




-1 ± 1 


K-TT-K+ 


(2.4567 ± 0.0003) x 10" 


-2 


(2.4572 ±0.0001) x 10" 


-2 


-0.05 ±0.03 




(2.2694 ± 0.0002) x 10" 


-2 


(2.2696 ±0.0001) x 10" 


-2 


-0.02 ±0.03 




(2.0877 ± 0.0001) X 10" 


-3 


(2.0878 ±0.0001) X 10" 


-3 


-0.001 ±0.002 




(2.5421 ± 0.0003) x 10" 


-2 


(2.5419 ±0.0001) X 10- 


-2 


0.03 ±0.03 




0.13077 ± 0.00001 




0.13078 ±0.00001 




-0.1 ±0.2 




0.13204 ± 0.00002 




0.13201 ±0.00001 




0.3 ±0.2 




(4.6175 ± 0.0009) x 10" 


-2 


(4.6160 ± 0.0008) X IQ- 


-2 


0.15 ±0.12 


CLEO model for Three Pions 


vr vr vr^ 


1.2605 ± 0.0001 




1.2604 ±0.0001 




1 ± 1 


^-^OyrO 


1.2702 ± 0.0001 




1.2702 ± 0.0001 




0± 1 


Two Pions and a Photon 


vr vr'^7 


(1.2708 ± 0.0001) X 10- 


-2 


(1.2707 ±0.0001) X 10- 


-2 


0.01 ±0.01 






4 


pions 






2vr-vr°vr+ 


0.95120 ± 0.00009 




0.95112 ±0.00006 




0.8 ± 1.1 


3vr°7r- 


0.24453 ± 0.0002 




0.24452 ± 0.0002 




0.1 ±2 






5 


pions 






vr-4vr° 


(3.6192 ± 0.0010) X 10- 


^2 


(3.6185 ± 0.0007) x IQ- 


-2 


0.07 ±0.12 


37r"27r+ 


(4.2745 ± 0.0010) x 10" 


-2 


(4.2746 ± 0.0006) x IQ- 


-2 


-0.01 ±0.12 


2vr"2vr°vr+ 


(1.1883 ±0.0004) X 10- 


-1 


(1.1884 ±0.0003) X 10- 


-1 


-0.1 ±0.5 



Table 2 



The partial widths for the three-, four- and five-meson decay modes of the r cal- 
culated using Herwig±+ and TAUOLA. In order to compare the results of the two 
programs, the masses of the decay products in TAUOLA were set to the Herwig±± 
values, as was the Cabibbo angle. For the five pion decays the p form factors were 
not included in the decay of the to. 

3. 6. 2 CLEO Model for Three Pions 

This is the implementation of the model of [31] for the weak current for three 
pions. This model includes p mesons in both the s- and p-wave, the scalar a 
resonance, the tensor /2 resonance and scalar /o(1370). The form factors for 
the vr°vr°vr~ mode are given in [31] and the others can be obtained by isospin 
rotation. 
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Fig. 9. Differential distribution for the mass of a) 7r^7r^7r+ and b) vr+vr" produced in 
the decay r~ tt~ tt^ tt^ Ur ■ The sohd hne is the model of |9|22|30j and the dashed 
line that of [31j . 





1.00 1.25 
m„V",r-/GeV 



0.75 1.00 



Fig. 10. Differential distribution for the mass of a) vr'^vr'^Tr" and b) tt^tt^ produced 
in the decay t~ tt^tt^tt^ i^r- The solid line is the model of [9|22|30j and the dashed 
line that of |31j . 

In this case the running width for the ai is calculated using the current for 
the TT^Ti^Ti^ and Tr^vr'^Tr" modes of the ai together with an s-wave KK* con- 
tribution. The running width is shown in Fig.fTSh for the default parameter 
values and our calculation is in good agreement with that in TAUOLA. 

The partial widths for the tt^tt^tt+z/^ and —>■ 71^71^71^ 1/^- are compared 

with those from the CLEO version of TAUOLA in Tabled The mass distribu- 
tions of the 7r~7r~7r"'" system and the vr+vr^ subsystem for the r~ ti'tt'tt^Uj. 
are shown in Fig.[9l The mass distributions of the hadronic system and the 
7r~7r° subsystem for the r~ — > Ti^n^Ti^Vr decay are shown in Fig. [10 There is 
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a) K TT b) TT 




1.2 1.4 1.6 1.8 0.6 O.a 1.0 1.2 



Fig. 11. Differential distribution for the mass of a) vr and b) vr produced 
in the decay r~ ^ i^~7r~Er^z^T-. The soHd hue is the model of [9|30j and the dashed 
line that of [25j. 



a) K°TT K° b) TT K 




1.2 1.4 1.6 1.8 0.6 0.8 1.0 1.2 



Fig. 12. Differential distribution for the mass of a) K^tt^ and b) tt^ produced 
in the decay r~ — > K^ir' K^Vt-. The solid line is the model of |9|30j and the dashed 
line that of 1251 . 



good agreement with TAUOLA for this model which in general gives a higher 
mass for the hadronic system and a slightly broader distribution for the masses 
of the subsystems which include a rho resonance than the model of [9ji22a30j . 
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Fig. 13. Differential distribution for the mass of a) K~tt^K^ and b) K^n^ produced 
in the decay r~ — > K~tt^ K^v-r- The sohd hne is the model of [9|30j and the dashed 
line that of 1251. 




Fig. 14. Differential distribution for the mass of a) K and b) tt^K produced 
in the decay r -n^iT^K^Ur- The solid line is the model of |9|30j and the dashed 
line that of [25]. 

3.6.3 Model for modes including Kaons 

Like the model of [30] the model of [25] is designed to reproduce the correct 
chiral limit for tau decays to three mesons. However, this model makes a 
different choice of the resonances to use away from this limit for the decays 
involving at least one kaon and in the treatment of the Ki resonances. 

The form factors for the different modes are given in [25]. The same form 
of the Oi Breit-Wigner is used as in Section 13.6.11 with the mass and width 
taken from [25]. The running width we calculate using this model with this 
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Fig. 15. Differential distribution for the mass of a) K^tt^tt'^ and b) vr+vr^ produced 
in the decay r~ — > K~7r~Tr^iyT-- The soHd hne is the model of [9|30j and the dashed 
line that of 1251. 




Fig. 16. Differential distribution for the mass of a) 7:~K^tt^ and b) tt~ produced 
in the decay r~ K^tt^i^t- The solid line is the model of |9|30j and the dashed 

line that of 1251 . 



choice of parameters is shown in Fig.[8)D and is in good agreement with the 
parameterisation given in fI2] . 

As well as the different choice for the resonances contributing to the various 
decay modes this model differs from that of [SUj in the treatment of the Ki res- 
onances. While the model of [30] assumes that only the ii'i(1400) contributes 
the model of [25] assumes that both the i^i(1270) and ^1(1400) contribute 
with different relative contributions to modes including i^*7r and Kp interme- 
diate states. All the parameters are taken from [25j . 
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Fig. 17. Differential distribution for the mass of a) vr tt^t] and b) vr vr'' produced in 

the decay t~ — > -K'TT^rjUr for the model of [9|30j . 

a) CLEO Model b) Four Pion Model 

2.0 I I I I I I I I I I I I I I I I I I I I I I I I I I I 0.5 I — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — I 




0-5 1.0 1.5 2.0 2.5 3.0 0.5 1.0 1.5 2.0 2.5 



Fig. 18. Running width of the oi in the models of [STj and [32]. a) shows the running 
width in the model of [31] with the dotted line showing the contribution of KK*, 
the dashed line showing the contribution of 7r~7r~7r''', the dot-dashed line showing 
the contribution of TT^ir^ir^ and the solid line the total running width, b) shows the 
running width in the model of |32j . 



The mass distributions for the total hadronic mass and the mass of the 7t~K~^ 
subsystem are shown in Fig. [11] for the decay r~ —>■ K~n^ K^Ur- The mass 
distributions for the total hadronic mass and the mass of the 7t~K^ system is 
shown in Fig.[T2]for the decay r~ — > K^n~ K^Ut-. The only major difference 
with the model of [30] for these decay modes is due to the different parameters 
for the ai and the inclusion of higher K* resonances. The mass distributions for 
the hadronic system and the K~ti^ subsystem for the decay r~ K~tt^K^Vj. 
are shown in Fig. [13] In this case the model of [25] includes K* resonances in 
the K~TT^ and tt^K^ subsystems which are not presented in the model of [30] 
giving different results for the masses of these systems and leading to a higher 
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a) K gTT Kg b) 




Fig. 19. Differential distribution for the mass of a) K^ir and b) K^n produced 
in the decay t~ — > K^'k~ K^Vr- The sohd hue is the model of [25] and the dashed 



line is the result of the model of |9j[30J . 




Fig. 20. Differential distribution for the mass of a) K^t:^ and b) K^ir^ produced 
in the decay r~ K^-k~ K^Vt- The solid line is the model of [53] and the dashed 
line is the result of the model of [9(30] . 

total hadronic mass for this decay mode. 



The mass distributions of the hadronic system and the tx^K~ subsystem in the 
decay r~ — > tt^'k^K~Vj. are shown in Fig.[T31 In this case the main difference 
between the models of [25] and (30] is due to the inclusion of the i^'i(1270) 
in the model of [25j which gives the two peak structure seen in the hadronic 
mass spectrum. The mass distributions for the hadronic system and the vr+vr" 
subsystem for the decay r~ — ^ i^~7r~7r are shown in Fig. [151 As with 
the Ti^TT^K' the inclusion of the i^i(1270) gives a different spectrum for the 
mass of the hadronic system when compared with that from the model of 
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Fig. 21. Differential distribution for the mass of a) vr n^^ and b) vr 7 produced in 
the decay t~ 'K~iT^^Vr for the model of [9]. 

|30] . The CLEO measurement of this mode [33] favours a significant i^i(1270) 
contribution but with a much larger width than the model of [25] for the Ki. 

The mass distributions for the hadronic system and the tt^K^ subsystem for 
the decay -k^ K^-k^Vt are shown in Fig.[T6l As with the K~tt^K^ the 

presence of K* resonances in the tt~K^ and K^ir^ systems give a different 
mass distribution for these systems, which together with the inclusion of the 
7^1(1270) gives the different result for the mass of the hadronic system. 

The mass distribution for the Kg7T~Kg and Kgir^ systems in the decay 
Kg-K^KgUj. are shown in Fig. [191 The mass distribution for the KgTT~K^ 
are K^tc^ systems in the decay r~ —>■ Kg7r~ K^Ur are shown in Fig.[20l The 
results from the model of [9f30] are also shown, where the and are 
allowed to decay with equal probability to Kg and without changing the 
hadronic current]^ The models are in reasonable agreement for the shapes of 
these distributions. 



3. 7 Two Pions and a Photon 



The branching ratio for the decay r~ — >■ ujtt'Vt- is 1.95% [34] • The majority 
of this decay is modelled as an intermediate state in the four pion current 
described below. However there is a 8.90% [31] branching ratio of the uj into 
7r°7 which must also be modelled. We do this using a current for 7r^7r'^7 via an 
intermediate uj. The hadronic current for this mode, together with the masses, 
widths and other parameters, are taken from [9]. 

^ The dashed lines in Figures [19] and [20] are therefore identical. 
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Fig. 23. Differential distribution for the mass of a) and b) produced 

in the decay t~ — > Tr'^Tr'^Tr'^Tr"!/.^ for the model of [32]. 



The partial width for the r~ — 7r'^7r~7 mode is given in Table. [2l The mass dis- 
tributions of the hadronic system and the n^'f subsystem, which contains the 
u resonance, are given in Fig.[2Tl There is good agreement between Herwig++ 
and TAUOLA for both the partial width and the shapes of the distributions. 
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Fig. 24. Differential distribution for the mass of a) vr vr'^ and b) vr 4it^ produced in 
the decay t~ A-K^iT~Vr for the model of [35j. 



a) TT^TT b) 37T 27T^ 




0.4 0.6 0.8 1.0 1.0 1.2 1.4 1.6 



Fig. 25. Differential distribution for the mass of a) vr^vr and b) Svr 27r+ produced 
in the decay t~ — > 37r~27r"'"j/.r for the model of |35j . 

c?. 5 Four Pions 



We use the model of |32j ^1 to model the decay of the r to four pions. The model 
is based on a fit to e~^e~ data from Novosibirsk. A Breit-Wigner distribution 
with a running width, Eqn. lA.2] is used for the a and u resonances with the 
running widths taken from [lOj. The more complicated form of [10] is used for 
the p which, apart from the choice to normalise to — 1 at = 0, is the same 
as the form of [23] given in Eqn. lA.4] 



It should be noted that there were a number of mistakes in this paper which were 
corrected in |10j. 
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Fig. 26. Differential distribution for the mass of a) vr+vr vr" and b) 2tt~-k+2-k^ pro- 
duced in the decay r~ — > A-K^iT~Vr- The sohd hne is the model of j35] neglecting 
the /9 propagators in the u decay and the dashed line includes these propagators, 
by default we do not include the p propagators. 

A Breit-Wigner distribution with a running width is used for the oi with the 
running width calculated in [32]. The Herwig-f-f calculation of the running 
width is compared with that from TAUOLA in Fig.fTSb. There is reasonable 
agreement between the two calculations apart from in the threshold region. 

The partial widths for the two modes are given in Table [2l The mass distri- 
butions for the hadronic system and the tt^tt+tt" subsystem, which contains 
the uj resonance, are shown in Fig.[22]for the decay -k^-k^-k^'k^Vt-. The 

mass distributions for the hadronic system and the tt^ti~ subsystem are shown 
in Fig. [23] for the decay r~ tt^'k^t:^'k~Vt-. There is good agreement between 
Herwig++ and TAUOLA for both the partial widths and mass distributions. 



3.9 Five pions 

We use the model of [35] which includes pu and pa intermediate states, via 
the ai meson to model the five pion decay modes of the r. 

The partial widths for these decay modes are compared with those from 
TAUOLA in Table[2l The mass distributions of the tt~7t^ subsystem and the 
hadronic system are shown in Fig. [21] for the decay Att^tt^Ut-. The mass 

distributions of the tt+tt" subsystem and the total hadronic mass for the decay 
T~ — * 37r^27r+z/T- are shown in Fig. [25] The mass distributions of the vr+Tr^vr" 
subsystem, which includes the ut resonance, and the total hadronic system are 
shown for the decay r~ — >■ 27r~27r°7r"'"z/.r are shown in Fig.[26] In all cases there 
is good agreement between the results of Herwig++ and TAUOLA for both the 



24 



partial widths and the shapes of the distributionsjZ] 



4 Modelling of Tau Decays 

In our approach the branching ratios for a given decay mode are specified as 
input parameters rather than calculated from the theoretical models used to 
give the distributions of the decay products, which gives us the ability to adjust 
the branching ratios to the experimentally-observed values. In general we have 
taken the decay modes and branching ratios for tau decays from [Mj. In [51] 
a set of basis modes is used, the branching ratios for which are constrained to 
sum to unity. In most cases we have used these modes. However, in some cases 
due to our modelling of the decays we have combined modes. We have also 
included some of the non-basis modes, which forces us to adjust the branching 
ratios slightly from those in [31] in order to ensure that the branching ratios 
still sum to one. For the higher multiplicity modes where charged hadrons 
are observed, but the different rates for kaons and pions are not known, we 
assume that the hadrons are pions. The decay modes, branching ratios and 
currents used are summarized in Table[3l In addition we have chosen to use 
the masses from [51] for the external particles, although we use the default 
choice of the models we are using, which often come from experimental fits, 
for the intermediate resonances. 

The branching ratios for the leptonic modes are taken from [51] and the dis- 
tribution of the decay products is modelled using the LeptonNeutrinoCurrent 
described in Section 13.31 The decays to a single charged meson are modelled 
using the ScalarMesonCurrent, described in Section [3TTI with the branching ra- 
tios taken from [34]- The Ktt modes are modelled using KPiCurrent, described 
in Section 13.51 and the parameters from the fit of [21] . The remaining two 
meson modes are modelled using the TwoMesonRhoKStarCurrent, described in 
Section 13. 4[ and the branching ratios from |j34j , the resonance parameters for 
the different modes are described in Section 

The three-pion modes are modelled using the ThreePionCLEOCurrent, as de- 
scribed in Section [3. 6. 2[ together with the branching ratios from [34J. We take 
the branching ratios for modes with three mesons, at least one of which is 
a kaon, from [51| and use the KaonThreeMesonCurrent, as described in Sec- 
tion I3.6.3[ to model the distribution of the decay products. 



^ As the u tt'^'k't:^ decay is below the threshold for /jvr production it can 
be modeled as either a contact interaction or via the pn intermediate state. The 
parameters of [35j are chosen to give the correct relative rates for the different 
intermediate states without the p propagators and we take this as our default choice. 
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Table 3 



Branching 


Decay Mode 


Current 


Ratio 






0.178345 
0.173545 




LeptonNeutrinoCurrent 
LeptonNeutrinoCurrent 


0.108924 
0.006885 


K-Ur 


ScalarMesonCurrent 
ScalarMesonCurrent 


0.254890 
0.008957 

0.004491 
0.001513 
0.000263 




TwoMesonRhoKStarCurrent 
KPiCurrent 

KPiCurrent 
TwoMesonRhoKStarCurrent 
TwoMesonRhoKStarCurrent 


0.092370 
0.089813 


TT^ TT^ IT' Vr 


ThreePionCLEOCurrent 
ThreePionCLEOCurrent 


0.003757 
0.003292 
0.000555 




KaonThrccMcsonCurrcnt 
KaonThreeMesonCurrent 
KaonThreeMesonCurrent 


0.001519 
0.001518 
0.001087 

0.000235 
0.000235 


KlK%TT-Ur 


KaonThreeMesonCurrent 
KaonThreeMesonCurrent 
KaonThreeMesonCurrent 

KaonThrccMcsonCurrcnt 
KaonThrccMcsonCurrcnt 


0.044435 
0.010313 
0.001762 


'K~ T:~^ T:~ TT^ Vr 


FourPionNovosibirskCurrent 
FourPionNovosibirskCTirrcnt 
TwoPionPhotonCurrcnt 


0.004935 
0.001744 
0.000957 
0.000834 


7r~ 7r~ TT"*" vr'^Tr'^ z^T- 

TT TT TT TT TT Z^'t 


FivePionCurrent 
ThreeMesonDefaultCurrent 
FivePionCurrent 
FivePionCurrent 


0.000225 

0.000145 
0.000135 
0.000118 


U>TT~ TT^TT^Ur 
UJTT~ TT~ TT^ Vt 


Phase Space 
Phase Space 
Phase Space 
Phase Space 


0.000400 
0.000397 
0.000307 
0.000280 
0.000238 
0.000225 




Phase Space 
Phase Space 
Phase Space 
Phase Space 
Phase Space 
Phase Space 


0.000297 
0.000059 


K-K+ TT-TT^Vr 


Phase Space 
Phase Space 



Decay modes and branching ratios used for r decays in Herwig++ together with 
the model of the hadronic current used for the decay. 
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We split the observed r~ uj-rVj- rate into its dominant pieces, i.e. uj 7r+7r~7r' 
and uj — > 71*^7. The three pion mode is included with the rest of the four pion 
tau decays and is modelled using the FourPionNovosibirskCurrent, described 
in Section [231 whereas the 7r°7 is modelled using the TwoPionPhotonCurrent, 
described in Section 13771 The branching ratios are taken from [31]. In general 
the treatment of non-dominant modes is often a problem when we include 
resonances as intermediate particles in the hadronic currents. However, in the 
other cases, such as the uj component to the five pion decays, the contribution 
of the other modes is smaller and comparable with other neglected tau decay 
modes and is not included. 

The five pion decays for the tau are modelled using the FivePion Current, de- 
scribed in Section [3^ As with the four-pion modes we include the uj contribu- 
tion with the rest of the five-pion modes as the current includes modelling of 
the intermediate uj contribution. In this case we neglect the sub-dominant uj 
modes due to their smaller contribution. The branching ratios are taken from 

m- 

The inclusion of the ri?)Tx and uj^iti modes is sufficient to saturate the observed 
six pion rates from [34| and we therefore include these modes using a phase- 
space distribution for the decay products. 

We also include a number of modes with one kaon and three pions, or two 
kaons and two pions and K*ri using a phase-space distribution for the decay 
products and the branching ratios from [51] . The largest branching ratio for 
one of these modes is Kuj with a branching ratio of 0.0004 and the sum of the 
branching ratios for these modes is less than three per mille. This means that 
they make a relatively small contribution and if necessary the simulation can 
be improved by implementing hadronic currents for these modes. 

As the decay modes, the decay models and branching ratios are specified in a 
data file these can easily be changed given new experimental results or better 
modelling of the hadronic currents. 



5 Other Applications 

The hadronic currents obtained from r decay can be used to describe other 
decays which occur via the weak current. It is an important test of the new 
structure of decays and weak currents in Herwig-|--|- that the currents can be 
used to simulate these decays. 

One obvious application is to use the currents, together with the naive factor- 
ization approximation, to simulate the weak hadronic decays of bottom and 
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a) Width 



b) Branching Ratios 
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Fig. 27. The a) width and b) branching ratios of the chargino as a function of the 
mass difference between the chargino and Hghtest neutrahno. The sohd hne shows 
the sum of the leptonic modes, i.e. xt Xi^^^^e and Xi ~^ XiM^^^m- "^^^ dashed 
hne shows the pion mode xt ~^ Xi'^^ and the dotted shows the xt ~^ Xi'^"''^^ 
mode. The dot-dashed hne shows the sum of the xt ~^ Xi^tt modes. The long 
dashed hne shows the sum of xt ~^ Xi{4, 5}7r, Xi" — *■ Xi'^'^'^^l^ and xt ~^ X?"^^ 
modes. The long dot-dashed line shows the sum of the modes involving kaons. 



charm hadronslZ] Here we will consider another possible use of these currents, 
the simulation of the weak decay of BSM particles where there is a small mass 
difference between two of the particles. One such example of this is the decay 
of the lighter chargino, which is almost mass degenerate with the lightest neu- 
tralino, in AMSB models. In most studies of these models only the leptonic 
and single-pion modes are included, in some more sophisticated studies the 
two- and three-pion modes were also considered [37f38] . 



Using the hadronic currents from tau decays it is easy to extend the simula- 
tion of BSM decays in [15] to automatically calculate the partial widths and 
branching ratios for decays with X —>■ YW^ where the is highly off-shell. 
These decays are calculated by using the Vertices for the BSM model, which 
encode the Feynman rules, and the hadronic current. The partial widths are 
normalised in order to give the correct branching ratios for the r with the 
modes which are modelled using a phase-space distribution neglected. This 
is a reasonable approximation given the relatively small contribution of these 
modes. As the partial widths depend on both the overall mass scale, as well 
as the mass difference, the branching ratios, while often close to, are not nec- 
essarily the same as those of the r when the mass difference is equal to the 
tau mass. In practice, this approach is normally only used for small values of 
the mass splitting and matched to the perturbative calculation with outgoing 
quarks at a mass difference below the charm threshold. 



We will use this approach in the simulation of hadronic decays in Herwig++ |36j . 
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We illustrate this by considering the decay Xi ~^ Xi^ at the AMSB pa- 
rameter point SPS9 [39]. The SUSY parameters were generated using SOFT- 
SUSY ^0|. The mass difference between the lightest neutralino and chargino 
is very sensitive to higher-order corrections, this can be seen by the large 
difference between the values obtained using SOFTSUSY (1.167GeV) and 
ISAJET7.58 [41j (0.167 GeV), which was used in [39j. Given the uncertainty 
in the mass difference, and to show dependence of the results upon it, we 
present the width and branching ratios as a function of the mass difference. 
Fig. [23 between the neutralino and chargino by keeping all the other SUSY 
parameters fixed to the values given by SOFTSUSY and varying the chargino 
mass. As can be seen, for larger values of the mass difference the multi-pion 
modes (> 3 pions) dominate, although at lower values they typically make a 
small contribution. 



6 Conclusion 

We have described the simulation of tau decays in the Herwig++ event gen- 
erator. A number of tests of the simulation have been performed giving us 
confidence in the implementation of the various hadronic currents and accu- 
racy of the simulation. 

This new simulation makes use of the factorization of the matrix element and 
the inheritance mechanism of C-|--|- to produce a simulation which is easy to 
extend by implementing new hadronic currents. The interface feature of the 
ThePEG framework [6], on which Herwig+-t- is based, allows easy access to all 
the parameters of the hadronic currents, allowing different experimental fits 
to be used without changing any code. 

The range of hadronic currents that are included together with the simulation 
of spin correlation effects in all tau decays gives a state-of-the-art description of 
tau decays. The new code structure will make this simulation easy to maintain 
and develop in the future. 
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A Breit-Wigner Distributions 



A number of different forms of the Breit-Wigner distribution are used in the 
various hadronic currents. The simplest form of the distribution, which is 
rarely used, is a fixed width Breit-Wigner distribution 

M^-iMT , 

= W—Etf- 

where s is the virtual mass squared of the resonance, M is the mass of the 
resonance and F is the width. As with many of the forms we use, this is 
designed so that B(0) = 1. 

In most cases we make use of a more sophisticated form including a running 
width for the particle, r(s): 

The details of this form of the Breit-Wigner distribution depend on the choice 
of the running width. If the decays of the resonance are dominated by a simple 
two-body decay, the running width is given by 




s p{M^)' 

TP- wave = r — f ) , (A.3b) 

(A.3c) 



where F is the physical width and p{s) is the momentum for the 2-body decay 
of a particle of virtual mass s in the rest frame of the particle. The angular 
momentum of the decay products is denoted as s-, p-, and d-wave for spin 0, 1, 
and 2 systems. For particles where the dominant decay is a three-body mode, 
for example the ai meson, a more complicated form of the running width is 
often used. In this case we still have BW(0) = 1, as the running width is zero 
below the threshold for the decay modes. 



In some cases for the p meson we use the form of Gounaris and Sakurai [23j|, 
which includes a running mass correction derived from an effective range for- 
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mula for the p-wave vrvr scattering phase, 



BW(^ 



M2 + dMT 



M^-s + T^ {A;2 [h{s) - /i(M2)] + kljh>{M^){M^ - s)} - iMT (^) ^ 

(A.4) 

where M is the mass of the p and F is its width. The function k{s) is defined 
to be 



k{s) 



-s — mi 



s > Ami 



(A.5a) 
(A.5b) 



where we have used /cm = k{M'^) and kg = k{s) in the definition of the Breit- 
Wigner distribution. The function h{s) is taken to be 



2 k 
I p cot 



Ami — s 



s > Am^ 
s < Ami, 



and h'(s) = The constant d 



IT k\.j 



2m^ 



2nk 



M 



mlM 
Txk^ 



M 



is defined to give BW(0) = 1. 
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